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I.   INTRODUCTION 

'  The  use  of  a  factorial  set  of  treatment  combinations  has 
now  become  a  widely  accepted  means  of  investigation.   The  prin- 
ciple of  factorial  experimental  designs  has  been  elaborated 
primarily  for  agricultural  field  experiments,  but  has  been 
found  valuable  in  many  other  types  of  experimentation.   Yates 
(1937)  has  described  many  designs  for  factors  at  two  and  three 
levels.   By  these  methods,  many  factors  may  be  introduced  into 
one  experiment,  but  the  total  number  of  plots  or  other  experi- 
mental units  required  for  a  comprehensive  test  is  large.   For 
example,  eight  factors  each  used  at  two  levels  require  2$6    plots 
for  one  replication,  and  six  factors  each  at  three  levels  re- 
quire 729  plots.   It  is  seldom  practicable  to  carry  out  experi- 
ments involving  so  many  plots.  A   further  device  for  reducing 
the  number  of  plots  is  that  of  fractional  replication,  whereby 
certain  treatment  interactions  are  assumed  negligible  and  only 
a  selection  of  all  possible  treatment  combinations  is  used  in 
the  experiment. 

In  this  report,  an  attempt  has  been  made  to  formulate  the 
general  principles  of  construction  of  fractional  asymmetrical 
designs  of  the  type  2m   and  3n  factorials  which  have  been  used 
to  obtain  the  factorial  replicated  designs  of  the  mixed  type 
2m  x  3n.   The  idea  behind  this  work  is  to  estimate  the  main 
effects  and  the  two-factor  interactions  (assuming  higher-factor 
interactions  to  be  absent)  with  as  feiv  assemblies  as  possible 
and  without  too  much  computation. 


II.   GENERAL  OUTLINE  OP  CONSTRUCTION  OF 

THE  FRACTIONALLY  REPLICATED 

EXPERIMENTAL  DESIGNS 


2.1  Treatment  Combinations  and  Their  Response 

Let  there  be  m  factors  A-^  A2,  .  .  •  >  Am  each  at  two 
levels  and  n  factors  B^,  B2,  .  .  . ,  B   each  at  three  levels. 
Then  there  are  2m  x  3n  treatment  combinations. 

The  symbol 

xl    x?  xm   ,  zl   ,  z?  zn 

al  >  a2  *  •  •  -,  %  ,  blX,  b2%  .  .  .,  bnn        (2.1.1) 

denotes  the  treatment  in  which  the  factors  A-,,  A2,  .  .  .,  A 
occur  at  the  levels  xlf    x2,  .  .  . ,  x^  x±  =   0,  1  ( i  =  1,  2, 
.  .  . ,  m)  and  the  factors  B-^  B2,  .  .  . ,  B   occur  at  the  levels 
z±,    z2,    .  .  .,  zn;  Zj.  =  0,  1,  2  ( j  =  l,  2,  .  .  .,  n)  . 

For  simplicity  the  treatment  combination  (2.1.1)  will  be 
called  an  assembly  and  denoted  by 

{x±,   x2,  .  .  .,  xm,  z±,    z2,  .  .  .,  zn)  (2.1.2) 

In  accordance  with  standard  convention,  denote  an  assembly 
and  the  mean  response  to  an  assembly  by  the  same  symbol.   Thus 
if  Y(x^,  x2,  .  .  .,  xm,  Zj,  z2,  .  .  .,  zn)  is  an  observed  re- 
sponse corresponding  to  an  assembly  (2.1.2),  then 
E[Y(x1,  x2,  .  .  .,  x^,  zl9    z2,  .  .  .,  zn)] 

=  [xlt    x2,  .  .  .,  xm,  z1,  z2,  .  .  .,  zn) 

where  E  stands  for  'expectation'. 
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2.2  Fractionally  keplicated  Asymmetrical  Factorial  Experi- 
mental Designs 


Let  S  be  a  column  vector  of  u  assemblies,  say, 

(XH>  xi2>  •  •  '»    xim) 
21*  "X:22,  *  *  •'  x2m 


S  = 


(*ul'  *u2> 


*um} 


which  also  is  an  array. of  some  strength  in  Euclidean  Geometry 
EG-(m,  2),  by  Barnard  (1936),  i.e.,  factorial  experimental  de- 
sign with  m  factors  each  at  2  levels.   Let  T  be  a  column  vector 
of  v  assemblies,  say 


T  = 


(zll>  z12>  * 
'z21*  z22'  * 


•'  zln} 

•'  z2n) 


(zvl>  zv2>  *  *  •'    zvn^ 


which  is  again  an  array  of  some  strength  in  EG(n,  3).   Then 
S  &  T  will  mean  the  symbolic  direct  product  of  two  column  vec- 
tors of  assemblies  S  and  T  and  is  taken  as 


k 


S  &  T  = 


(x-q,  x12, 


'  Xlm>  zll'  z12' 


^\l'  *u2>  *  '  •'  xum'  zll'  z12' 


ln; 


^xll*  x12>  *  *  •>    xlm>  zvl>  zv2>  •  '  *>  zvn' 


'  zlnJ 


(-^ul*  xu2»  •  •  •*    xum>  zvl>    zv2'  •  •  •>    zvn) 


which  is  a  (uv  x  1)  column  vector  of  mixed  assemblies.   It  can 
also  be  looked  upon  as  a  fractionally  replicated  design  con- 
sisting of  uv  assemblies  of  a  2m   x  3n  design.   (Ref . :  Mac- 
Duff  ee  (12)) 


^m 


2.3  Effect  of  the  2m  Factorial  Experiment 


-.m 


Any   interaction   of   a    2      experiment   will   be   denoted   by 


•m 


12  m 


(k±  =  0,    1;    i   =   1,    2,    .    .    .,    m) 


(2.3.1) 


It  also  includes  the  main  effects  and  the  average  response  of 
all  assemblies.   If  k-  =  1  and  the  rest  of  the  k!s  (i  =  2.  3. 
.    .    .,  m)  are  zero,  (2.3.1)  represents  the  main  effect  A-,.   If 
k±  =   k2  =1  and  the  rest  of  the  k^s  (i  =  3,  4,  .  .  . ,  m)  are 
zero,  then  the  (2.3.1)  represents  the  two-factor  interaction 
A-jA2>  and  similarly  for  other  factors. 


Let 

1 

-l" 

"l 

o(O)"" 

H  = 

= 

1 

1 

1 

c(l)_ 

(2.3.2) 


then 


H 


-1 


'  1/2    1/2" 
-1/2    1/2 


where  the  constants  c(x),  (x  =  0,  1) ,  are  defined  by  (2.3.2) 
Let 


U^  =TJ&U&U&...&TJ 


denote  the  product  of  U  taken  m  times.   Then  we  shall  define 
all  the  interactions  in  (2.3.1)  by  the  matrix  identity 


1  ' 

& 

1 

_al_ 

_a2_ 

am 


=  EM 


i 
A™ 

_  m_ 
(2.3.3) 


or 


I  ' 

a 

i 

LAi_ 

W 

& 


with  the  convention 


Q 


m 


=  H" 


(m) 


1  " 

& 

1    " 

K81- 

_a2_ 

&.  .  .& 


l 

am 


1-1  =  1,   I  Ak  =  Ak  =  Ak  I,   A°  =  I 

The  expression  for  the  typical  effect  A^1  A^2    .  .  .  A.km 


m 


is 


Akl  Ak2 
Al  A2 


.  .   A%*  =  J)7T    (-1  +  3i)xi  (1  +  a^1-^: 
where  the  divisor  D  =  (1/21),  i  =  1,  2,  .  .  .,  m. 
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For  example,  if  m  =  3, 

I  =  1/8  (1  +  ai)(l  +  a2)(l  +  a3) 

A-l  =  1/8  (-1  +  ax)(l  +  a2)(l  +  33) 

A1A2A3  =  1/8  (-1  +  ai)(-l  +  a2)(-l  +  a3) 

where  I  is  the  average  of  the  mean  responses  of  all  assem- 
blies.  (Ref . :  Bose  and  Conrior  (5).) 

Assuming  three  or  more  factor  interactions  as  negligible, 
(2.3.3)  reduces  to 

,  xl   x2        x™        3- 
(al   a2   •  •  ■  am  )  =1+^7  c(xi)Ai 

i=l 


+  Z-   c(xi)  c(xi,)  AiAi, 


i<i 


where  c(x±) ,    x±   =   0,  1,  i  =  1,  2,  .  .  .,  m  are  as  given  in 

(2.3.3)  and  2-       means  the  summation  over  all  pairs  of  indices 
i<i' 

i,  i'  for  i  <  i<  (i,  i'  =  1,  2,  .  .  .,  m). 
2.I4.  Effects  of  the  3n  Factorial  Experiment 

Any  interaction  of  a  3n  experiment  will  be  denoted  by 

_WT     WO  Wr, 

V  V  •  '  •  V-  (2.^.1) 

(w,  =  0,  1,  2;  j  =  1,  2,  .  .  .,  n). 


;s  a 


When  only  one  w  is  different  from  zero,  (2.1^.1)  represent: 
main  effect.   It  is  a  linear  effect  if  w  =  1  and  a  quadratic 
effect  if  w  =  2.   When  two  w's  are  different  from  zero,  (2.1^.1) 


represents  a  two-factor  interaction.   It  is  a  linear  x  linear, 
linear  x  quadratic  effect  according  as  the  couplet  of  non-zero 
w'a  is  (1,1),  (1,2)  or  (2,1)  and  (2,2). 
Let 


K  = 


1 
1 
1 


-1 

0 
I 


1 
-2 

1 


1  dx(0)  d2(0) 
1  d-^1)  d2(l) 
1   d1(2)    d2(2)j 


(2.1J..2) 


Then 


K 


-1 


1/3 

1/3 

1/3 

-1/2 

0 

1/2 

1/6 

-1/3 

1/6 

where  the  constants  dw(z);  z  =  0,  1,  2;  w  =  1,  2  are  defined 
by  (2.2J..2).   Then  we  shall  define  all  interactions  in  (2.1^.1) 
by  the  matrix  identity 


1  " 

~i~ 

bl 

& 

b2 

k 

jt 

or 


I 

I  " 

Bl 

& 

B2 

LBi_ 

_BL 

=  K<*)  . 


a 


=  K 


(-n) 


I 

i ' 

Bl 

& 

B2 

LBL 

B2 

1      ' 

1     " 

^1 

& 

bg 

-bij 

31 

(2.1J..3) 


a 


with  the  convention 


1-1=1,   IB 


w 


jW  _  DW 


Bw  =  BWI,   B' 
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W-i    W  ;  Wn 

The  expression  for  the  typical  effect  B-,x  B..  '  .    .  .  B  u  is 
Bjl  Bf2  .  .  .  s^n 
=  D7T  (l+bi+b2)f0(zj)(-i+b2)fl(2j)(1„2b  b2)f2(zj) 

•  J     J  J  J     J 

J 

where  D=(21x3°x6i2)~,n  =  the  number  of  B  factors  with 
exponent  w.,  nQ  +  n-^   +  n2  =  n,  (j  =  1,  2,  .    .    . ,  n)  ; 

f0(Zj.)  =  l/2(zj  -  l)(zj  -  2),  ^(Zj)  =  Zj(-Zj  +■  2),  and 
f2(Zj)  =  1/2  Zj(Zj  -  1). 

For  example,  if  n  =  2, 

I  =  1/9(1  +  bx  +  bi)(l  +  b2  +  b2) 

B^g  =  1/12 (-1  +  bf)(l  -  2b2  +  b|). 

As  before,  I  is  the  average  mean  response,  and  the  effects  de- 
fined here  do  not  obey  the  convention  that  the  sum  of  the  posi- 
tive  coefficients  is  unity.   B-jB2  is  the  interaction  of  the 
linear  component  of  B1  and  the  quadratic  component  of  B2. 
(Ref . :  Bose  and  Conner  (5).) 

Assuming  three  or  more  factor  interactions  as  negligible,  - 
(2.ij..3)  reduces  to 

W=l  j=l       J     J 


>       I  dw(Zj)  dw,(z  )  B^, 


where  dw(z-),   w  =  ±,  2,      z.  =  0,  1,  2,   j  =  1,  2,    .  .  .,    n, 
are  given  in  (2.1|.2),  and  £_   means  the  summation  over  all 

pairs  of  indices  j,  j'  for  j  <  j'  (j,  j'  =  1,    2,    .  .  . ,  n) , 

and  2l   means  the  summation  over  all  pairs  of  indices 
w,w? 

w,  w1  (w,  wT  =  1,  2) . 


2.5  Effects  of  the  2m   x  3  Factorial  Experiment 


Any  interaction  in  this  experiment  will  be  denoted  by 
,kl  Ak2 


km  w^_   W2 
A-i   Ao   .  •  .  A    B-i   Bp   • 


B 


_w 

n 


n 


(2.5.D 


k-i   kp 
which  may  be  regarded  as  the  symbolic  product  of  Aj  Ap 


vm 


W]_   W2 


w 


n 


.  .  .  Amm  and  B]_"L  B2   .  .  .  Bn  .   If  kj_  =  1  and  w,-  =  1  and  the 
rest  of  (k,  w) ' s  are  zero,  (2.5.1)  represents  the  linear  effect 
A^_  Bj  of  the  ith  factor  of  the  2m   factorial  and  the  jth  factor 
of  the  3n  factorial.   If  k*  =  1  and  w-  =  2  and  the  rest  of 

J 

(k,  w)'s  are  zero,  then  (2.5.1)  represents  the  quadratic  effect 
All  interactions  in  (2.5.1)  will  be  defined  by  the  matrix 


A.jB   of  the  same  factors  (i  =  1,    2,    .  .  . ,  m;  j  =  1,  2,  .  .  . ,  n) 


identity 


1 
al 


1 
a2 


1 

1  ' 

1  " 

a 

bl 

& 

bp 

a^ 

_  m 

A 

A 

& 
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=  Eim)   a  K(n>  . 


a. .  .a 


I 

■  1 

i 

_Am 

a 

Bi 

_Bi 

a 

b2 

R2 

_  2_ 

a.,  .a 


A/ 


& 


~I  ' 

i 

I  " 

S 

Bl 

a 

B2 

m 

A 

R2 
_  2_ 

a 


i 

a 

Bn 

B2 

n 

-  tt-^)  .  Tr-(n) 


=  E 


K 


~1  " 

a 

1  " 

_al_ 

_a2_ 

a.,  .a 


i 

l  ' 

l 

_v 

a 

bi 

a 

b2 

^2 
_bl_ 

A 

a. .  .a 


(Ref.:  Bose  and  Connor  (5).) 


The  effect  of  A^l  A^  .  .  .  AJ%  B™1  B™2  .  .  .  Bwn  is 
dTT  [(-l+ai)Xi(l+a1)1"Xi(l+bj+b^)fo(ZJ^-l+bj)fl(2J) 


(l-2b,+b^)f2(zj) 


J  "J 

ni+m    n. 


_i  * 


,  n-i+m    nn    no  -1 

where  DM21   x  3  °  x  6  2)    L,    ^   =  the  number  of  B  factors 

with  exponent  Wj>  nQ  +  nx  +  n2  =  n,  ( J  =  1,  2,    .    .  . ,  n) ,  and 

fo(zj^  fl^zj)»  f2(zj^  are  defined  below  (2.2j..3).   (Ref.: 
Connor,  ¥.  S.  (7).)   If  m  =  3,  n  =  2 

A^Bf  =  1/96  (-l+a1)(l+a2)(l+a2)(-l+b^)(l-2b2+b2) 

where  I  denotes  the  average  of  the  mean  responses  of  all  mixed 
assemblies . 

Assuming  three  or  more  factor  interactions  as-  negligible, 
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(a-,1  a^'  .  .  .  amri  •  b-^  b2<  .  .  .  bnn) 

m  ^-  n 


m  ^~  2   n 

=  I   +    %_     c(xi)Ai  +  2.   c(xi)c(xi,)AiAi  +  2T   X  dw(xi)Bi 
i=l  i<i '  w=l  i=l 


+  2   X   dJ^K^z.,) 


2 
w  w ' 


W    I  '  w 
1   i-Ci  1         J 


J    J  J    w=l  i,j  J     J 

(2.5.2) 

where  c(x- ) ,  d  (z-)>  >£-  >   .2  ,    are  as  already  defined  earlier, 

and  5_  means  the  summation  over  all  pairs  of  indices  i,  j 

(  i  =  1 ,  2,  .  .  . ,  m ;  j  =  1 ,  2 ,  .    .    . ,  n). 

2.6  Fractional  Replicate  of  the  2m   Factorial  Experimental 
Design 

The  treatment  combination  of  the  2m  experiment  may  be 
represented  as  points  of  an  n-dimensional  lattice  with  axes 
x-|_,  x2,  .  .  • ,  xm  (as  given  in  section  2.1).   Each  axis  of  the 
lattice  will  have  two  points,  these  points  being  the  elements 
of  the  Galois  Field  GF(2).   Thus,  for  example,  the  treatment 
combinations  of  the  2   factorial  experiment  can  be  represented 
by  the  points  (0,  0),  (l,  0),  and  (1,  1).   The  two  factors  may 
be  denoted  by  Ax  and  A2.   The  points  (0,  0)  and  (1,  1)  both 
satisfy  the  equation  x-j_  +  x2  =  0,  and  comprise  a  1/2  replicate 
of  the  2   experiment  defined  by  the  identity  relationship 
I  =  A-j_A2.   The  symbol  A-|_A2  can  be  used  to  denote  the  set  of 
treatment  combinations  for  which  x,  +  Xp  =  0 .   Similarly,  the 
symbol  A  A   denotes  the  set  of  treatment  combinations  for 
which  x-,  +  x2  =  1. 
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Consider  m  factors,  A-^  A2,  .  •  •  >  Am  each  with  two  levels. 
Then  A-^1  A2  Ay  .  .  .  A^1  denotes  the  set  of  treatment  combina- 
tions for  which  ^xl   +   k2x2  +  •  •  •  +  ^m*™  =   i>  where  i  is  an 

element  of  the  Galois  Field,  GF(2).   A  subscript  i  can  be 
associated  with  every  effect  or  interaction  of  the  identity  re- 
lationship.  For  example,  the  treatment  combinations  of  the  four 
possible  l/ij.  replicates  of  the  2?    experiment  defined  by  the 
identity  relationship, 

i  =  A*i  42  A*3  „  A*i  A*k  Ap  =  42  43  A*k  k\s 

are  those  treatment  combinations  which  satisfy 

(1)  k^x-]_  +  k2x2  +  koXo  =  0 

klxl  +  *V%  +  k^x5  =  0 

k2x2  +  k^x^  +  k^  +  k^x^  =  0 

(2)  k1x1  +  k2x2  +  knXo  =  0 

klxl  +  \*\+  +   k5x5  =  1 

k2x2  +  k^x3  +  k^  +  k^X£  =  1 

(3)  kixi  +  k2x2  +  knx^  =1 

klxl  +  k4xi4.  +  k5^5  =  0 

k2x2  +  k3x3  +  k^x^_  +  k^x^  =  1 


(i|.)  k-^  +  k2x2  +  kjx^  =  1 
klxl  +  \^  +  k^x5  =  1 
k2x2  +  k3x3  +  k^x^  +  k^x^  =  0  , 


respectively, 
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Since  the  interaction  ApAoAj^ Ah  is  the  generalized  inter- 
action of  A^A2A^  and  AjAlA?,    the  subscripts  associated  with 
A2A3AkA5  in  the  four  fractional  replicates  may  be  obtained  as 
the  sum  of  the  corresponding  subscripts  of  A-jApAo  and  A-,AlAh. 

Thus  a  l/2r  fraction  of  the  2m  factorial  consists  of  the 
assemblies  (x^  x2,  .  .  .,  x^  satisfying  the  r  linearly  inde- 
pendent equations. 

Ls  =  aslxl  +  as2x2  +  •  •  •  +  asmxm  =  ds  (2.6.1) 

(ds  =  0,  1),  (s  =  1,  2,  .  .  .,  r), 

in  Galois  Field  (2),  where 

(asl-  as2>  •     -     •>    asm)  1   (0,  0,  •  •  -,0). 

Because  each  of  the  ds  can  assume  one  of  two  values,  0  or  1, 
there  are  2r  different  fractions  which  can  be  produced  in  this 
way.   Each  fraction  contains  2rn-r  treatment  combinations,  and 
because  the  different  fractions  do  not  have  any  treatment  com- 
binations in  common,  the  fractions  collectively  contain  all  2m 
treatment  combinations. 

We  shall  say  that  the  interaction, 

41  42    ■    ■    ■    A>  (2.6.2) 


corresponds  to  the  linear  form  k^  +  k2x2  +  .  .  .  +  k  x 
in  GF(2) . 

The  interaction  (2.6.2)  is  estimated  by  the  contrast 
such  that 


Ik 


1/2™  [(klXl  r  k2X2  +.;.+  k^  =  1) 

-  (klXl  +  k2x2  +  .  .  .  +  k^  =  0)]  (2.6.3) 

where  the  parenthesis  (  )  means  the  sum  of  the  observed  re- 
sponses of  assemblies  satisfying  the  equation  within  it.   If 
only  one  k  differs  from  zero,  (2.6.2)  represents  a  main  effect. 
If  two  k's  differ  from  zero,  it  represents  a  two-factor  inter- 
action, and  so  on. 

Consider  any  linear  form  L  which  is  not  of  the  form, 

klLl  +  k2L2  +  ■  •  •  +  Vr  (2.6.^) 

ka  =  0,  1;  s  =  1,  2,    .  .  .  r 

(kx,  k2,  .  .  .  kp)  ?£   (0,  0,  .  .  .,0) 

the  interactions  corresponding  to  the  linear  form, 
L  +  (k1L1  +  k2L2  +  .  .  .  +  kpLr) 

are  said  to  be  aliases  of  the  interaction  corresponding  to  L. 

It  is  known  from  the  theory  of  fractional  replication  that 
each  interaction  not  corresponding  to  any  linear  form  k-.L,    + 
k2L2  +  .  .  .  +  kpLr  is  a  member  of  one  and  only  one  alias  set. 

It  is  clear  that  the  assemblies  of  the  fraction  defined  by' ' 
(2.6.1)  which  satisfy  L  =  0  also  satisfy  L  +  (k1L1  +  k2L2  + 
.  .  .  +  krLr)  =  ^d!  +  k2d2  +  .  .  .  +  krdr  and  those  satisfying 
L  =  1,  satisfy  L  +  (k^  +  k2L2  +  .  .  .  +  kpLr)  =  1  +  k^ 

+  k2d2  +  .  .  .  +  krdr,  where     ksds  =  0  or  1.   The  rank  of 

s=l 
the  matrix  formed  by  the  column  vectors  of  coefficients  belong- 
ing to  effects  which  correspond  to  L  +  (k^-,  +  k2L2  +  .  .  . 
+  krLr)  is  1. 
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The  identity  relationships  for  the  fraction  defined  by 
(2.6.1)  are 

I  =  Gx  =  G2  =  .  .  .  =  Gp  =  G1G2  =  .  .  .  =  Gr.1Gr 

=  .  .  .  =  G-L&2  .  .  .  Gr  (2.6.5) 

where  the  G's  are  the  interactions  corresponding  to  the  linear 
form  (2.6.4) .  . 

The  relationship  (2.6.5)  is  useful  in  determining  the 
set  of  effects  aliased  with  a  given  effects. 

2.7  Fractional  Replicate  of  the  3n  Factorial  Experimental 
Design 

In  a  3n  design,  suppose  each  treatment  combination  to  be 
represented  in  the  conventional  fashion  by  a  code  containing  n 
digits.   The  Kth  digit  represents  the  level  of  the  Kth  factor, 
and  take  the  values  0,  1,  2  representing  the  three  levels,  as 
given  in  section  2.1^.   Then  each  such  code  may  be  thought  of  as 
the  coordinates  of  a  point  in  an  n-dimensional  factor  space, 
and  the  whole  design  consists  of  3n  such  points  arranged  in  a 
hypercubic  lattice.   This  geometrical  way  of  considering  the 
fractional  replicate  of  the  3n  factorial  design  agrees  with 
the  method  given  in  section  2.6. 

Thus  a  l/3q  fractional  replicate  of  a  3U  factorial  consists 
of  the  assemblies  (z^  z2,  .  .  .,  zj  satisfying  the  q  linearly 
independent  equations 

Mt  =  btlzx  +  bt2z2  +  .  .  .  +  btriZn  =  et  (2>?al) 

(et  =  0,  1,  2),  (t  =  1,  2,     .  .  .,  q) 
in  GF(3). 
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We  shall  say  that  the  interaction 


B™1   B£>  .  .  .  B^ 


corresponds  to  the  linear  form  w^z^  +  w2z2  +  .  .  .  +  w  z   which 
we  may  assign  to  the  two  orthogonal  components,  linear  and 

quadratic,  denoted  by  LCB™1  B22  .  .  .  B™n)  and  QiB™1   B22  .  .  . 

Bnn)  ,  respectively.   The  linear  effect  L(B11  B22  .  .  .  B™n)    is 
estimated  by  the  contrast 

1     r 

2  x  3n-l  L(W1Z1  +  W2Z2  +  •  •  •  +  wnzn  =  2) 

-  (w^Zj  +  w2z2  +  .  .  .  +  wnzn  =  0)J 
and  the  quadratic  effect  QtB™1  B™2  •  .  .  B™n) ,  by  the  contrast 

1     r 

(W-.Z-,  +  w9z9  +  .  .  .  +  w  z  =  0) 

2  x  3n  '- 

-  2(w1z1  +  w2z2  +  .  .  .  +  wnz  =  1) 
+  (wlZl  +  .  .  .  +  wnzn  =  2)] 

where  (  )  has  the  same  meaning  as  in  (2.6.3). 

Let  M  be  any  linear  form  which  is  not  of  the  form, 

w^  +  w2M2  +  .  .  .  +  wqMq  (2.7.2) 

wt  =  0,  1,  2;  t  =  1,  2,  .  .  .,  q; 

(wx,  w2,  .  .  .,  Wq)  ^  (0,  0,  .  .  .,  0). 

The  interactions  corresponding  to  the  linear  form  M  +  (w  M 

+  w2M2  +  .  .  .  +  wqMq)  are  said  to  be  aliases  of  the  interaction 
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corresponding  to  the  linear  form  11. 

It  is  clear  that  the  assemblies  of  the  fraction  defined 
by  (2.7.1)  satisfying  M  =  0,  1,  2  also  satisfy 

M  +  (w^  +  w2M2  +  .  .  .  +  w  M  ) 
r3-        S  q 

-  X   wtet*  E,  wtet  +  i.   H  wtet  +  2 
t^i  t=i  t=i 

respectively,  where  e.  =  0, '1,  2. 

The  identity  relationships  for  fraction  defined  by  (2.7.1) 
are 

11  =Gi  =   Q2=    •  •  ■  «*i"0{<4=  •  •  ■  =^.!G' 

...  =  ...=  G{G2Gj  .  .  .  G^2  ■  (2.7.3) 

where  the  (G')'s  are  the  interactions  corresponding  to  the 
linear  form  (2.7.2) . 

The  relationship  (2.7-3)  is  useful  in  determining  the  set 
of  effects  aliased  with  a  given  effect. 

2.8  Fractional  Replicate  of  the  2m  x  3n  Factorial  Experimental 
Design  r 

W.  S.  Connor  (7)  has  published,  an  example  of  his  general 
method  of  fractionating  2™  x  3n  designs.   The  example  is  a  half-' 
replicate  of  a  2^   x  32  design.   Briefly,  the  method  by  which  it 
is  constructed  is  as  follows.   The  2^    part  of  the  design  is 
split  into  two  conventional  half-replicates  S±   and  Sp.   The  32 
part  of  the  design  is  split  into  three  conventional  one-third 
replicates  T±,    T2,  and  T3 .   The  final  design  can  then  be  repre- 
sented  symbolically  by 
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S1T1  +  S2T2  +  S2T3 

where  S1T1,  for  example,  represents  all  possible  combinations 
of  treatment  combinations  in  Sx  with  treatment  combinations  on 
T1.   This  gives  a  design  of  36  out  of  a  total  of  72  trials.   For 
constructing  plans  for  the  2m   x  3n  arrangements  are  adjoining 
fractions  of  the  3n  factorial  to  fractions  of  the  2m   factorial 
in  such  a  way  as  to  preserve  two-factor  interaction  estimates. 
To  illustrate  the  method  of  construction,  consider  the  2^  x  32 
arrangement.   A  l/k   replicate  of  the  2^  x  32  can  be  constructed 
by  adjoining  the  combinations  of  a  1/3  replicate  of  the  32 
factorial  to  a  3/k-   replicate  of  the  2^  factorial. 

An  interaction  A^  A*2  .  .  .  A£n  B"l  B^2  .  .  .  B™n  may  be 

said  to  belong  to  two  linear  forms  k-jX,  +  k2x2  +  .  .  .  +  k  x 
in  GP(2)  and  w^  +   w2z2  +  .  .  .  +  w^  in  GP(3). 

Mi?  42  .  .  .  A^  b£  B^  .  .  .  ajn," 

and    Q(A^1  A^2  .  .  .  A"%  B»l  B|2  .  .  .  B£n,  . 

Let  the  equations 

kxxx  +  k2x2  +  .  .  .  +  k^  =  d 
wlZl  +  w2z2  +  .  .  .  +  WnZn  =  e 

be  satisfied  by  points  U±t    x,,  .  .  .,  ^  z±f    z,,  .  .  .,  ^ 
such  that.(Xl,  x2,  •  •  .,  xm)  satisfy  k^  +  k2x2  +  .  .  . 

+  ^  =  d  in  GF(2>>  and  <*1,  ^2,  .  .  .,  zn)  satisfy  wlZl  +  w2z2 
+  .  .  .  +  wnzn  =  e  in  GP(3),  then  (2  x  3n_1  x  2m) 


19 


• 

LU^'42  •   • 

Akm  Rwl  Rw2 
•    ^Sn     bl     b2 

contrast 

/klxl  +   •    • 

•    +   ^^  =   -1-] 

>']_Z]_    +     .     . 

•    +   wnzn  =   2' 

.    B   lL)    is   estimated  by   trie 


klxl   +    •    •    -    +   VSn  : 

W]_Z]_    +     . 


.    +  wnzn  =  0 


kixi  +  •  -  -  +  Vm 

wlzl  +    •    • 


knXn      +      . 


=   0, 


+   "m^  "   u\     _   /Klxl   +    '     '     '    +   "nrSi 
+  wnzn  =   2'         \W]_Z]_   +    .     .     .    +  wnzn 

and(2x3nx   2m)    x  Q(A^   A^2    .     .     .    A^  B?1  B^2    .     .    .    B™n) 


is  estimated  by  the  contrast 


"  /klxl  +  . 

•    •    +   ^m^m  =   1\ 

o( 

k-j_x-j_   +    .    . 

•    +   K*m   =   1) 

\WjZ-j_    +     . 

.     .    +   wnzn  =   0/           W-lZj   +    .     . 

•    +  wnzn  =   l' 

+  /klxl  +  • 

•    '    +  wnzn  =   2  _ 

- 

/klxl  +    • 

\W-j_Z-j_    +     . 

•     '    +   VSi  =  ° 
.    .    +  wnzn  =  ° 

^x-l  +  .  .  .  +  kmzm  =  0^    /k1x1   ■   .  .  .  +  kmxm  =  0 


f*lxl  +  •  •  •  +  Vm  =  Ox 
^•WjZ-l  +  .  .  .  +  wnzn  =  l' 


W-j_Z-|_  +  . 


+  wnzn  =  2 


where  (  )  has  the  same  meaning  as  in  (2.6.3). 

The  identity  relationships  for  the  fractional  replicate 
of  2m   x  3n  design  can  be  obtained  from  those  of  2m  and  3n  de- 
sign.  If  the  fraction  is  given  by  the  mixed  assemblies  ob- 
tained by  combining  the  assemblies  of  the  fraction  defined  by 
(2.6.1)  and  the  assemblies  of  the  fraction  defined  by  (2.7.1) 
as  in  section  2.2,  then  the  identity  relationships  are 
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I  =  Gn  =  G^ 


G1G2 


-   G±G2 


■    Gp  —  G^  -  Gp  — 

=  GiGJ  =  G2G2  = 

i  t  , 

GpG^Gp   .  •  .  Qq 


GiG2 

•  =  GrGq 


where  the  interactions  are  obtained  by  taking  all  products  of 
interactions  from  the  identity  relationships  for  the  fractions 
defined  by  (2.6.1)  and  (2.7.1). 

To  make  these  points  clear,  consider  a  fraction  of  a 
2C  x   3^  design  consisting  of  the  assemblies  (x-]_,  x2,  zlt    z2) 
obtained  by  taking  the  symbolic  direct  product  of  assemblies  of 
an  array  given  by  Xl  +  x2  =  0  in  EG(m,  2)  and  those  of  an  array 
given  by  z±  +   z2  =   0  in  EG(n,  3). 

Writing 

and  T- 


Si  = 


(0,  0) 

(1,  1) 


(1,  2) 

(2,  1) 


the  required  assemblies  are  given  by  the  definition  in  section 
(2.2),  thus 


S]_  &  T±   = 


(o,  0,  0,  0)' 

(0,  0,  1,  2) 

(0,  0,  2,  1) 

(1.  i,  o,  0) 

(1,  1,  1,  2) 

_(1,  1,  2,  1) 


(2.8.1) 


Let 
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P'    =|jl,    A±,    A2,    A±A2    ,    L(BX),    Q(B1),    L(B2)  ,    Q(B2),    LCB^g), 
Q(B-lB2),    KB-lB^),    Q(Bxb|),    L(A1B1),    qCAjlB!),    LfA^), 
Q(A-lB2),    KAgBi),    Q(A2BX),    L(A2B2)  ,    Q(A2B2)~j      . 


Then 


E   = 


"y(0,    0,  0,    0)' 

7(0,    0,  1,    2) 

y(o,  o,  2,  i) 

yd,  i,  o,  o) 

yd,  i,  i,  2) 

yd,  i,  2,  i) 


i  - 

t  _ 
J. 

i  - 

i 
i 

i 

•  p 


i 
i 


(2.8.2) 
1-1      1-1      1-1      1   -1" 
0      2-1-1      0      2-1-1 


1   -1   +1    -1      1-1      1-1  1   -1 

1-1+1      0-2      1      1-1  1      1 

1-1+1      1      1      0-2-1  1      0-2-1-1      0      2-1-1      0      2 

1      1+1-1      1-1      1-1  1-1      i-i      i-i      i    _!      x    _x      l 

11+10-211-11110-2110-211 
1      1+1      1      1      0-2-1      1      0-2      1      1      0-2      1      1      0-2 


The  identity  relationship  for  the  fraction  (2.8.1)  and  the 
other  aliased  sets  of  effects  are 


I  =  A3A2  =  B1B2 


[Al>  A 2] 

[3,  B2,  BlB|] 


(2.8.3) 


\_AjB-l,    AgB-L,  A^,  A2B2J 
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assuming  higher  factor  interaction  to  be  absent. 

Thus  the  six  estimable  functions  of  effects  are  one  each 
for  the  first  two  and  two  each  for  the  remaining  two  sets  of 
effects  in  (2.8.3) • 

From  (2.8.2),  the  normal  equations  estimating  them  are 


6 

-6 

-6 

6 

6 

6     6 

-6 

6 

6     6 

-6 

6 

-6   -6 

6 

6 

6" 

h-A  - 

Al 

6 

6 

A2 

A1A2 
Q(BXB2) 


Y(AX) 
Y(A2) 


"KD 
Y(A^2) 

Y-KB^) 
Y.QCB-LBg) 


1|  0  2-6  2-6 

0  12  -6  -6  -6  -6 

2  -6  ij.  0  Ij.  0 

-6  -6  0  12  0  12 

2  -6  ij.  0  l±     0 

•6  -6  0  12  0  12 


"l^b-l) 

~Y-L(B1) 

QCBi) 

Y.Q(BX) 

L?B2) 

Y-L(B2) 

Q(B2) 

Y-Q(B2) 

A          ? 

L.(B1B2) 

y.l(b1b|) 

A         2 

_Q(B1B2) 

y.q(b1b|) 

(1) 


(2) 


(3) 
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k  0  2   -6  Ij.     0  2  -6' 

0  12  -6   -6  0   12  -6   -6 

2  -6  ij.     0  2   -6  ^     0 

-6  -6  0   12  -6   -6  0   12 

4  0  2  -6  Ij.     0  2  -6 

0  12  -6  -6  o  12  -6  -6 

2  -6  ij.     0  2   -6  1±     0 

■6  -6  0   12  -6   -6  0    12 


QU^) 
L(A1B2) 
Q(A-,B2) 
L(A2BX) 
Q(A2B1) 
L(A2B2) 
Q(A2B2) 


J   L 


YUA^)" 

YQ(A1B1) 

YL(A1B2) 

YQU-^) 

YL(A2B!) 

YQ(A2Bi) 

YL(A2B2) 

YQ(A2B2) 

(4) 


(2.8.4) 


which  implies  the  estimability  of 


I  -  AXA2  -  L(BXB2)  +  QCB-LBg)  , 


(2.8.5) 


Ax  +  A2  , 

i|L(B1)  +  2L(B2)  -  6Q(B2)  +  2L(BXB2)  -  SQ^Bg)  , 

12Q(B1)  -  6L(B2)  -  6Q(B2)  -  6L(B1B2)  "  Sc^B*)  , 

^(A^)  +  2L(A1B2)  -  6Q(AXB2)  +  liL(A2B1)    +  2L(A2B2)  -  6q(A2B2)  , 

12Q(A1B1)  -  6L(AXB2)  -  6Q(AXB2)  +  12Q(A2B1)  -  6L(A2B2) 

-  6Q(A2B2)  . 

Thus  the  six  assemblies  of  22  x  32  factorial  design,  taken  as  '' 
above,  estimate  the  six  linear  functions  of  effects  given  in 
(2.8.5). 

The  Y  functions  on  the  right-hand  side  of  the  equations  in 

(2.8.4)  are  the  usual  C ' j  components  in  the  normal  equations 

A 


where 


C  C  P  =  C'j 

Sz)  =  c  P  . 
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III.   UNIVARIATE  LINEAR  MODEL  WITH 
FIXED  EFFECTS 

Let  the  fractional  design  consist  of  N  assemblies 
(N  -  2   x  3n)  with  y1}    y2,  .  .  . ,  yN  as  their  responses.   These 
N  responses  are  observable  quantities  with 

yr  =   0r  +  er     (r  =  1,  2,  .  .  .,  N)  (3.1) 

where  6r  is  the  true  value  of  yr  and  er  Is  an  error  in  observa- 
tion.  The  true  values  are  assumed  to  be  linear  functions  of 
main  effects  and  two-factor  interactions  (including  the  grand 
average  I  =  p.)  n,  |xx,  .  .  . ,  tx ( p-1)  ln  some  order,  which  are 
unknown  and  fixed  with  known  non- stochastic  coefficients  1, 
Crl>  Cr2-'  •  •  •»  Cr(p-1)  as  Siven  in  (2.5.2).   For  instance, 

.    Qr  =  V-  +   Crl^l  +  Cr2.u2  +  •  •  •  +  Cr(p-l)^(p-l)     (3.2) 
(r  =  1,  2,  .  .  .,  N) . 

Written  out  in  full,  the  coefficients  C    fr  =  1   2         n- 
v  =  1,  2,  .  .  ,t    (p-l)J  are  as  indicated  below. 

1.  If  jxv  corresponds  to  the  main  effect  A-  (i  =  1,  2, 

.  .  . ,  m)  of  the  2m  factorial  design,  (3.3) 

■1  if  yr  has  level  0  of  the  ith  factor 
1  if  jr   has  level  1  of  the  ith  factor. 


crv 


2.  If  nv  corresponds  to  the  two  main  effects  B  .,  B2. 

«J       J 

(j  =  1,  2,  .  .  .,  n)  of  the  3n  factorial  design, 
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th 


r-l   1  if  yr  has  level  0  of  the  j  n  factor 


rv 


0  -2  if  yr  has  level  1  of  the  jtn  fact 


or 


th 


■TV 


^  1   1  if  yr  has  level  2  of  the  jrn  factor. 

3.  If  [i-v  corresponds  to  the  interaction  Ai  Ai   (i,  /  i2 
=  1,  2,  .  .  . ,  m)  of  the  2'  factorial  design, 

(   1  if  yr  has  levels  0,  0  of  i-^  and  i2  factors 

1-1  if  yr  has  levels  0,  1  of  i-^  and  i2  factors 

Crv  =  \ 

\-l  if  yr  has  levels  1,  0  of  i-^  and  i2  factors 

\^1  lf  Jv   has  levels  1,  1  of  i1  and  i2  factors. 

2 

k-    If  M-v  corresponds  to  the  interactions  B.  B,  ,  B-  B- 

Jl  J2'   Jl  J2* 

BJlBJ2'  Bj1Bj2  (J'l  *   h   =   1>    2,  .  .  .,  n  of  the  3n 

factorial  design, 

/l  -1  -1   1  if  yr  has  levels  0,  0  of,  jx  and  j2  factors 

0   0   2  -2  if  yp  has  levels  1,  0  of  jx  and  j2  factors 

-1   1-1   1  if  yr  has  levels  2,  0  of  j±   and  j2  factors 

0   2  0  -2  if  yr  has  levels  0,  1  of  j±    and  j2  factors 

0   0   0   ii  if  yp  has  levels  1,  1  of  j±   and  j2  factors 

0-2   0-2  if  yr  has  levels  2,  1  of  ^  and  j2  factors 

-1-1   1   1  if  yr  has  levels  0,  2  of  ^  and  j2  factors  • 

0   0-2-2  if  yr  has  levels  1,  2  of  jx  and  j2  factors 

yi   1   1   1  if  yr  has  levels  2,  2  of  jx  and  j2  factors. 

5.  If  Hv  corresponds  to  the  interactions  AjB-,  A-B- 

(i  =  1,  2,  .  .  . ,  m;    j  =  1,  2,     .    .    . ,  n)  of  the  2m   x  3n 
design, 
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rv 


f  1  -1  if  yr  has  levels  0,  0  of  i  anu  j  factors 

0  2  if  yr  has  levels  0,  1  of  1  and  j  factors 

-1  -1  if  yr  has  levels  0,  2  of  i  and  j  factors 

-1  1  if  yr  has  levels  1,  0  of  i  and  j  factors 


0  -2  if  y  has  levels  1,  1  of  i  and  j  facto 


rs 


1  1  if  y  has  levels  1,  2  of  i  and  j  factors 


Thus  the  model  (2.5-2)  can  be  written  as 


Y  = 


£(Nxl)  ~  °(Nxp)  "  U(pxl)  +  l(Nxl) 


where   C  =  ((1,  Cpv) )  as  defined  in  (3.3), 

M'  =  [y->    V-if    \*-2>    •  •    •>    ^(p-1)]  > 
E(e)  =  0 


and 


V(e) 


o-h 


N 


( <T~     unknown) 
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IV.   MULTIVARIATE  LINEAR  MODEL  WITH 
FIXED  EFFECTS 


In  a  fractional-  design  with  N  assemblies  out  of  2m  x  3n 
(N  <  2   x  3  ),  there  are  N  experimental  units  to  be  observed. 
The  problem  is  univariate  or  multivariate  according  as  we  make 
one  observation  or  more  than  one  observation  on  each  experi- 
mental unit. 

Suppose  we  make  "q"  observations  ( q  >  1)  on  each  experi- 
mental unit,  then  the  model  analogous  to  that  in  (3.1)  is 


Yr  =  1 


y.   +  crlii1  +  cr2ii2  +  ...  +  op(p.1)ii(p.1)  +  er  (i+.l) 


where  each  j_p,  jx,  £±,    .    .    .,   jx^   ^  ,  ep  is  a  (q  x  1)  column 
vector,  ja,  ]j,v  (v  =  1,  2,  .  .  . ,  p-1)  are  assumed  to  be  fixed 
and  unknown,  and  e  is  NID  ( 0( q  x  1) ,    (  q  x  q) ) .   The  matrix  of 
constants  G  =  ((1,  crv) ) ,  r  =  1,  2,  .  .  . ,  N;  v  =  1,  2,  .  .  ., 
(p-1)  is  the  same  as  in  the  univariate  case. 

The  model  (if- 1)  can  be  compactly  written  in  the  form 


Y»(N  x  q)  =  C(N  x  p)  .  ja(pxq)  +  e«( 


Nxq), 


(4.2) 


where 


(Nxq) 


12 


In. 
q 


i 

l 


*2 


(Nxq) 


_-N_ 


1 

1 
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^(pxq) 


Hi 


1 
1 


-   (1) 


(q) 
^1 


q 

We  shall  assume  (i)  l(qx^)  to  be  an  unknown  positive  defi- 
nite matrix,  (ii)  N  >  p  +  q  where  p  is  the  rank  of  C. 
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V.   ESTIMATION  OP  EFFECTS  ( LEAST-SQUARE  METHOD) 

From  the  model  (3. It),  we  know  the  sum  of  squares  due  to 
error  is 

e'e  =  (Y  -  Cy.)  '  (Y  -  Cjx)  . 

The  value  of  jx  that  minimizes  e'e  is  given  by  the  solution  tc 

9    .  ■       ' 

(e'e)  =0. 

d  M 

We  get 

*  A 

(e'e)  =  2CY  -  2  C'Cjj,  =  0  . 

The  least-square  estimate  of  jx  , 

£  =  (C'C)"1  C'Y  , 

is  the  best  linear  unbiased  estimate  of  ix,  and  we  get 
var.  (jj)  =  (C'C)"1^2  . 

The  unbiased  estimate  of  <j~2   based  on  the  least-square 
estimate  of  jj,  is  given  by 

A     Y'(IN  -  C(C'C)"1  C')Y 
<T      =   _ 

N  -  p 

All  the  p  effects  represented  by  the  effect  vector  jx  are 
estimable  only  if  rank  C  =  rank  C'C  =  p;  otherwise  some  of 
them  are  non-estimable.   (Ref.:  Graybill  (10).) 
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VI.   TESTING  OP  SIGNIFICANCE  AMD 
CONFIDENCE  BOUNDS 

We  again  assumed  that  J3/N  -,%  is  NID  (0(Nx-i  )  >  <J~   ^m)  and  on 

the  basis  of  information  supplied  by  the  N  responses  y-i,  y?, 

.  .  . ,  y^  certain  test  procedures  about  the  nature  of  effects 
are  given. 

6.1  Single  Hypothesis  (Student's  t-test) 

p 

Consider    H  •  u.  =  0 

O   -C-0 

H:  u  ^  0 


where  [i   denotes  a  certain  specified  effect.   Let  /jl  be  the 
least-square  estimate  of  y.     with  var.  ([j,  )  =  d.O~    ,    where  d  is 
some  constant  and 


2«(iN  -  c(c»c)_1  c)i 

N  -  p 
(Ref . :  Graybill  (10) .) 


=  cr 


=  cr 


Then  an  *<  level  test  of  the  above  hypothesis  is  given  by 


A  (i.e.,  reject  H  )  if  t  = 

0(j1,  72>    •••*  yN)  =  \  'Vd 

^0  (i.e.,  accept  H  )  otherwise 


1*0 


>  t 


where  t^  is  such  that 


e[^(7!,  y2,  ...,  yN)/H0]  = 
6.2  Multiple  Hypothesis  (F-test) 


(0  <  -<-<!). 


Le 


t    (ji*)  '  =  (^  n2,  .  .  .,  nr) 


be  the  row  vector  of  r  effects  (r  <   p)  with 
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A.«.x  ,  /A         A 


t£ 


,^r) 


Then  an  -<   level   test   for   the   hypothesis 


0 

1 

0 

0       .     . 

.     .     .       0 

0       . 

,     .       0 

0 

0 

1 

0       ,     . 

.     .     .       0 

0      . 

.     .      0 

0 

0 

0 

1    .   . 

.     .     .      0 

0       . 

.       0 

0      0      0      0      .     . 
1  r 


0      10      ...      0_ 
P   -    (r+1) 


'   ^(pxl) 


or 


"  ^(rxl  *  -(rxl)    , 

HQ  =   (GH  =  £*)    =  0 
Ha   =   (Gy.  -  %*)    j*  0 


is   given  by 


Tl    (i.e.,    reject  HQ)    if  P  >  F^" 

#(yi,  i2>  •••'  y^  =  \ 

\0    (i.e.,    accept  H   )    elsewhere    , 


where 


P  = 


(N   -    p)    y_>    LC^C'C)"lc'    "    G(G'G)    LC-Jj_ 
r  y_r[lN   -   C(C'C)-1C  JZ 


and  F^    is   such  that 


Ef^yi,  y2,    •    .    .,  yN)/H0_    = 
(Ref . :    Graybill    (10) .) 


(0  <  «■<  1) 
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6.3  Confidence  Bounds  for  Estimate  (Single  Hypothesis) 


Given  an  estimable  effect  y.Q    (say),  the  probability  state- 


ment 


M^^o 
V  dO"2 


$** 


*(. 


based  on  student's  "t"  distribution  can  be  inverted  to  read 

rrf%  -   t,  jJ~&  s  ,oS  %o  +  t,V^}  -  i  -  <. 

6.I4.  Confidence  Bounds  for  Estimate  (Multiple  Hypothesis) 

The, variance  ratio  statistic,  F,  in  (6.2)  can  be  used  to 
give  a  probability  statement 

-  p    (£*  -  ji*)  '  [gCC'C)"^']"   (2-::-  -  H#) 


N 


z'[i  -  C(C'C)-1C'] 


<  p 


=  !-=<, 


which  implie 

P 


f[cjU)'   [ocoio-^i]"1  (£*)] V2 

F+  '  (z'  [i  -  c(cc)-1  cjv)~'1//2 

11/2 
1/2 


N   -    p 


(it*)'    [g(C'C)-1   G']"1    (M-::-) 
(£*)'     [G(C-C)-1   G'l"1    (£*) 

*77-^  *  (I't1  -  c(c.c)-1c]zJ 


1/2 


>   1   -   *    . 


(Ref.:    Roy,    S.    N.    and  Bose,    R.    C,    ( 1>+) ,    Roy   S.    N.    {!$)  . 
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The  inequality  in  the  parenthesis  provides  a  confidence  interval 
for  the  parametric  function 


i  -1 


(u*)<  [G(C'C)-1  G'j    (j±*) 


ll/2 


with  a  confidence  coefficient  ^  1  -  ^. 


3k 


VII ;   SUMMARY 

The  method  of  construction  is,  in  brief,  as  follows.   We 
obtain  t  fractions  of  the  2   complete  factorial  out  of  r  frac- 
tions of  the  3n  complete  factorial.   We  then  associate  each 
fraction  from  the  2   with  a  fraction  from  the  3n.   The  associa- 
tion is  such  that  every  treatment  combination  in  the  fraction 
from  the  2   is  adjoined  to  every  treatment  combination  in  the 
fraction  from  the  3n,    thereby  forming  treatment  combinations 
in  the  fraction  from  the  2m   x  3n. 

Section  2  contains  the  standard  method  for  obtaining  frac- 
tional factorial  designs  from  the  2m,  3n,    and  2m  x  3n  complete 
factorial.   An  example  is  given  to  describe  the  method  of  con- 
struction of  the  mixed  type  of  2m   x  3n  fractionally  replicated 
factorial  design.   In  sections  3  and  l+,    there  are  univariate 
and  multivariate  linear  model  analyses.   Section  f?  contains  a 
description  of  least  squares  estimation.   The  last  section  is 
to  test  hypotheses  that  are  whether  a  specified  (or  each  of  a 
set  of)  effects  is  significantly  different  from  zero,  or  the 
effects  within  a  given  set  are  all  equal. 
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The  use  of  a  factorial  set  of  treatment  combination  has 
now  become  a  widely  accepted  means  of  investigation.   The  prin- 
ciples of  factorial  experimental  design  have  been  elaborated 
primarily  for  agricultural  field  experiments,  but  have  been 
found  valuable  in  many  other  types  of  experimentation.   Yates 
(1937)  has  described  many  designs  for  factors  at  two  and  three 
levels.   By  these  methods,  many  factors  may  be  introduced  into 
one  experiment,  but  the  total  number  of  plots  or  other  experi- 
mental units  required  for  a  comprehensive  test  is  large.   For 
example,  eight  factors  each  used  at  two  levels  require  256  plots 
for  one  replication,  and  six  factors  each  at  three  levels  re- 
quire 729  plots;  it  is  seldom  practicable  to  carry  out  experi- 
ments involving  so  many  plots.   A  further  device  for  reducing 
the  number  of  plots  is  that  of  fractional  replication,  whereby 
certain  treatment  interactions  are  assumed  negligible  and  only 
a  selection  of  all  possible  treatment  combinations  is  used  in 
the  experiment. 

In  section  2,  an  attempt  has  been  made  to  formulate  the 
general  principles  of  construction  of  fractional  asymmetrical 
designs  of  the  type  2m  and  3n  factorials  which  have  been  used 
to  obtain  the  factorial  replicated  designs  of  the  mixed  type    ' 
2m  x  3n.   To  illustrate  the  method  of  construction,  consider  a 
half-replicate  of  a  23  x  32  deslgn<   The  ^    papt  Qf  ^  ^.^ 

is  split  into  two  conventional  half-replicates  Sx  and  S2.   The 
3   part  of  the  design  is  split  into  three  conventional  one-third 
replicates  Tx,  T2,  and  T3 .   The  final  design  can  then  be  repre- 
sented symbolically  by  S-^  +  S2T2  +  S^ .   Where  S^,  for 


example,  represent?  all  possible  combinations  of  treatment  com- 
binations in  S]_  with  treatment  combinations  on  T-,.   This  gives 
a  design  of  36  out  of  a  total  of  72  trials.   In  section  3  and  l\., 
there  are  univariate  and  multivariate  linear  model  analysis. 
Section  5  contains  a  description  of  least  square  estimation. 
The  last  section  is  to  test  hypotheses  that  are  whether  a 
specified  (or  each  of  a  set  of)  effects  is  significantly  dif- 
ferent  from  zero,  or  the  effects  within  a  given  set  are  all 
equal.   The  idea  behind  this  work  is  to  estimate  the  main  effects 
and  the  two-f actor  interactions  (assuming  higher  factor  inter- 
actions to  be  absent)  with  as  few  assemblies  as  possible  and 
without  too  much  computation. 


